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^ Abstract 

We report the observation of the superconducting proximity effect in nanoribbons of a candidate 
topological insulator (Bi2Se3) which is interfaced with superconducting (tungsten) contacts. We 
observe a supercurrent and multiple Andreev reflections for channel lengths that are much longer 
than the inelastic and diffusive thermal lengths deduced from normal state transport. This suggests 
^ that the proximity effect couples preferentially to a ballistic surface transport channel, even in the 

presence of a coexisting diffusive bulk channel. When a magnetic field is applied perpendicular to 
^ the plane of the nanoribbon, we observe magnetoresistance oscillations that are periodic in magnetic 

field. Quantitative comparison with a model of vortex blockade relates the occurrence of these 
I> oscillations to the formation of Pearl vortices in the region of proximity induced superconductivity. 
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I. INTRODUCTION 



Proximity-induced superconductivity in superconductor-normal conductor (SN) junctions 
has a long history of theoretical and experimental studj^. Recently, the superconducting 
proximity effect has attracted renewed theoretical attention within the context of topological 
insulators (TIs), materials wherein topologically protected, spin-polarized surface states are 
created by the combination of strong spin-orbit coupling and time reversal symmetry^ ^. In- 
terfacing a TI with a conventional superconductor is of fundamental interest for a variety of 
reasons. At such interfaces, theory predicts the formation of zero-energy mode quasiparticles 
that are condensed matter analogs of elementary fermionic excitations envisioned by Majo- 
rana but not yet observed in Naturd^. Furthermore, the "locking" of spin and momentum 
in TI surface states raises important theoretical questions about the nature of the proximity 
effect at Tl/super conductor junctiond^ and recent experiments suggest that the measure- 
ment geometry could have a significant influence on the observed phenomenaP. Evidence 
for the proximity effect was provided in an early study^ that interfaced a superconductor 
with Bii_^.Sb2:, a material now recognizecP^'^ as a 3D TI. More recentlj^, a supercurrent 
was observed in thin exfoliated samples of another TI (Bi2Se3). Bulk superconductivity has 
also been seen in compounds derived from a parent Tp3IIll_ However, systematic studies of 
proximity-induced superconductivity in candidate TIs are still in their infancy. 

In this paper, we discuss measurements of the proximity effect in mesoscopic Bi2Se3 chan- 
nels interfaced with superconducting tungsten (W) leads. Our principal aim is to report the 
observation of two new experimental results. First, we show that even in non-ideal TI sam- 
ples with bulk conduction, ballistic transport in the surface states can manifest through the 
persistence of a supercurrent and multiple Andreev reflections over significantly longer dis- 
tances than the phase breaking and diffusive thermal lengths deduced from (bulk dominated) 
normal state transport. Second, at temperatures above the onset of complete superconduc- 
tivity, we observe magnetoresistance (MR) oscillations that cannot be understood using 
conventional scenarios such as fluxoid quantization or the Aharonov-Bohm effect. Instead, 
the data are quantitatively explained using a recent modeP^ that relates MR oscillations 
in superconducting channels to the "Weber blockade" of Pearl vortices. The combined 
observation of proximity induced superconductivity and vortices in a Tl-superconductor 
configuration could be relevant in the ongoing search for Majorana fermiond^Ml^ 
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II. SYNTHESIS AND CHARACTERIZATION OF BI2SE3 NANORIBBONS 



The Bi2Se3 nanoribbons studied here were synthesized via gold catalyzed vapor-liquid- 
solid mechanism using a horizontal tube furnac^^*^. The source material, Bi2Se3 shots 
(Alfa Aesar, 99.999%), was placed in the center of the hot furnace where the temperature 
was kept at ~ 530 °C during the growth. A Si (100) substrate dispersed with gold catalyst 
particles (~ 20 nm in diameter) was placed downstream at the cold region zone. The tube 
was flushed with Ar several times before the growth to remove residual oxygen. A 60 seem 
Ar flow was kept at 1 Torr at the base pressure of ~ 10 mTorr during the growth. The 
furnace was cooled down to room temperature after a 1 hour and 30 minute growth. Figure 
la shows a scanning electron microscopy (SEM) image of an as-grown sample. A typical 
growth usually produces nanoribbons with thickness 60 < t < 100 nm, width 200 ^ w < 500 
nm and length 2 < L < 30 fim. Figure lb shows a high resolution transmission electron 
microscopy (TEM) image from the edge of a typical nanoribbon. The interplanar distance 
along the nanoribbon growth direction is 0.21 nm, which is consistent with the Bi2Se3 lattice 
constant [a = b = 0.4140 nm, c = 2.8636 nm). Figure Ic is a selected area electron diffraction 
pattern from the same ribbon: all the nanoribbons we have studied so far show a growth 
direction along [1120]. The clear lattice fringes and the electron diffraction pattern indicate 
that these ribbons are single crystals with little disorder. 

We obtained additional confirmation about the crystalline phase of the samples using 
Raman spectroscopy of individual Bi2Se3 nanoribbons. Figure Id shows a typical room 
temperature Raman spectrum from a single Bi2Se3 nanoribbon supported over one of the 
holes in a TEM grid; the data were taken with 514.5 nm excitation and the geometry ensures 
that the back scattered light only originates from the nanoribbon of interest. The incident 
radiation was polarized parallel to the nanoribbon growth direction and the scattered ra- 
diation was unpolarized. The measurements were performed with a very low incident laser 
power (~ 10 /iW with a spot size of ~ 1 /im) to avoid sample overheating. Three character- 
istic Raman modes are observed: ~ 72 cm~^, out of plane vibration mode; ~132 cm~^, in 
plane vibration mode; and ~ 173 cm~^, out of plane vibration mode. All these modes are 
consistent with those observed from bulk Bi2SeaP'. 

We fabricated electrical transport devices by transferring as-grown Bi2Se3 nanoribbons to 
a Si substrate with a 1 /im thick Si3N4 insulating layer and then depositing superconducting 
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tungsten and/or normal platinum (Pt) electrodes using a dual beam focused ion beam (FIB) 
etching and deposition system. The tungsten contacts are contaminated with carbon and 
gallium, and consequently have a high superconducting transition temperature (4 K< Tc< 
5 K) and a large critical field {He > 70 kOe)ff^. The correspondingly large superconducting 
energy gap (A ~ 0.7 meV) allows good energy resolution in differential conductance mea- 
surements. Furthermore, the large critical field lends such devices to explorations of possible 
vortex formation in regions of proximity-induced superconductivity. The FIB method also 
has the advantage of forming good ohmic contacts with Bi2Se3 nanoribbons, albeit at the 
price of localized damage in the contact area. TEM measurements of test structures (Fig. 
2a) show that the nanoribbons remain in the single crystal Bi2Se3 phase after FIB depo- 
sition, apart from the region immediately underneath the tungsten contacts. The typical 
electrode is about 200 — 400 nm wide and ~ 50 nm high. X-ray energy dispersive profile 
scans in scanning TEM studies (Fig. 2b) shows that the tungsten spreads about 250 nm 
from the visible edge of the contact. We conservatively estimate that the contacts have a 
resistance-area product RA ~ 10^^'^ Qm^. To rule out the possible shorting between closely 
spaced tungsten contacts, we also carried out a control measurement by depositing two tung- 
sten strips with a visible edge-to-edge separation of 460 nm on a Si3N4 substrate; we found 
that the (two-probe) resistance of this configuration was larger than 0.5 MQ at T = 500 
mK, well below the superconducting transition temperature of the tungsten contacts. To 
further rule out insidious spreading effects that could occur in a nanostructure geometry, we 
also examined (undoped) ZnSe nanowires with tungsten contacts and again found no signs 
of superconducting shorts below the transition temperature of the contacts. 

We measured a total of 6 Bi2Se3 nanoribbon devices with tungsten contacts; their char- 
acteristics are summarized in Table 1. In addition, we also measured two devices with 
non-superconducting (Pt) contacts. Although the nanoribbons were all fabricated in the 
same synthesis run, their dimensions and electrical properties (resistivity) show an inho- 
mogeneous distribution. We are unable to determine the carrier density and carrier mobility 
because the nanoribbons are too narrow to carry out Hall effect measurements; however, 
the electrical resistivity of all the devices (~ 1 mfi-cm) is of the same order as typical values 
observed in thin films of Bi2Se3, suggesting similar carrier densities (n ~ 10^^ cm~^). We 
measured the bias-dependent differential conductance in single Bi2Se3 nanoribbon devices 
using a dc current source and a lock-in amplifier at a frequency of 97 Hz over a temperature 
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range 0.5 K < T < 6 K and in magnetic fields up to 80 kOe. The measurements were per- 
formed in botli four probe and two probe geometries; tlie latter measurements are carried 
out in a "pseudo-four-probe" scheme, with two wires on each contact pad. 

III. TRANSPORT MEASUREMENTS WITH NORMAL METAL (PT) CON- 
TACTS 

Before we discuss transport measurements of Bi2Se3 nanoribbons with superconducting 
electrodes, we first study transport properties of two devices with non-superconducting 
electrodes in the four probe geometry. Figures 3a and 3b plot the temperature dependent 
conductivity for device G and H respectively, showing a ln(T) variation, consistent with 
2D quantum corrections due to weak anti-localization and electron-electron interaction^. 
Detailed fitting of cr vs T becomes more complicated because of the many fitting parameters: 
as we show below, fitting the magnetoconductivity to the weak localization theory suggests 
that there are multiple channels participating in the quantum corrections. Since the slope of 
cr vs. ln(T) is strongly affected by the dependence of L^{T), it is not appropriate to assume 
all channels contribute equally, as this assumes the coherence length for each of the channels 
is the same as the others. Figures 3c and 3d show the magnetoconductivity for device G and 
H at T = 3.0 K (squares) and T = 2.0 K (circles) in fields perpendicular to the nanoribbon 
plane. We fit the magnetoconductivity using the Hikami-Larkin-NagaokEpH equation: 

- A.w.(0) (* 4^/) - (i^/) ) • m 

where a depends on the dominant scattering mechanism and \l/(x) is the digamma function, 
a = 1 in the limit of weak spin orbit and magnetic scattering; a = —1/2 in the limit 
of strong spin orbit scattering and weak magnetic scattering; and a = when magnetic 
scattering is strong. For completeness, we also include the 2D correction from electron- 
electron interactions^^ given by: 

AaEEi(T, H) - AaEEi(r, 0) = -^F^g^iT, H), (2) 
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and is a function of the average of the static screened Coulomb interaction over the Fermi 
surface. The sohd hues in Figs. 3c and 3d are fits to the above equations, using as fitting 
parameters the phase breaking length L^, the parameter a and the screening parameter F^. 
For device G, the fits yield ~ 90.77 nm, a = —1.30 and F„ = 0.0 at T = 3.0 K and 
~ 137.2 nm, a = -1.34 and F^ = 0.0 at T = 2.0 K. For device H, the fits yield ~ 35.5 
nm, a = -1.36 and F^ = 0.32 at T = 3.0 K; ~ 47.0 nm, a = -1.29 and F^ = 0.36 at 
T = 2.0 K. Note that the values of the phase breaking length are somewhat smaller than 
typically seen in thin film samples of Bi2Se2r^. 

IV. TRANSPORT MEASUREMENTS WITH SUPERCONDUCTING (TUNG- 
STEN) CONTACTS: OBSERVATION OF A SUPERCURRENT AND MULTIPLE 
ANDREEV REFLECTIONS 

Figure 4a shows the temperature dependence of the two probe zero bias differential resis- 
tance {dV/ dl) of device A with width w = 600 nm and thickness t ~ 60 nm. The resistance 
is measured between contacts with edge-to-edge separation L = 1.08 |Um(SEM image in the 
inset). Spreading of the tungsten contacts implies a conservatively estimated Bi2Se3 chan- 
nel length of ~ 580 nm. The figure shows the onset of the proximity effect at T ~ 4.7 K 
when the W contacts become superconducting, eventually transitioning to a zero resistance 
supercurrent state at T < 2 K. Figure 4b plots the I — V characteristics of this device at 
different temperatures, showing a critical current Ic = 1.1 /uA at T = 500 mK. The product 
IcRn (where Rn is the normal state resistance) has a small value (~ 165 /iV << ^), simi- 
lar to observations in exfoliated Bi2Se3 layer^. As expected, Ic decreases with increasing 
temperature until the supercurrent is completely suppressed at T > 2 K. 

Figure 4c plots dl/dV vs. Va.tT = 500 mK, revealing "subharmonic gap structure" due 
to multiple Andreev reflections at V = where A is the gap of the W electrodes and n 
is an integei'^S^. The arrows identify a consistent series of subharmonic peaks in dl/dV 
corresponding to n = 2, 4, 8 (Fig. Id) for a superconducting energy gap A = 0.669 meV. 
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Using the BCS relation A = I.TG/cbTc, this corresponds to Tc= 4.41 K, which is close to 
Tc= 4.7 K obtained from the temperature dependent zero bias differential resistance. The 
absence of some values of n is not understood, but is not an uncommon occurrenc^^. The 
peak assignment in Fig. 4c results in an anomalous location of the prominent n = 1 peak 
at 1.194 meV instead of the expected location V = ^ = 1.338 meV (indicated by the arrow 
in Fig. 4c). We attribute this to heating at higher bias, noting that similar anomalies have 
been seen in other SNS devices^. The origin of the peak identified as A' = 0.932 meV is 
unclear; it could be a signature of a quasiparticle minigap Ag that is expected to arise when 
< (i.e. when the Thouless energy roughly equals the superconducting 
gap of the S contacts)ES. The data in Figs, la-d provide a clear indication of proximity- 
induced superconductivity in Bi2Se3. In Figs. 4e and 4f, we map out the detailed variation 
of dl/dV as a function of temperature and magnetic field, respectively. These plots show 
that - as expected from the temperature and field dependence of the superconducting gap 
of the tungsten electrodes - the multiple Andreev reflection features and Ic are quenched 
with increasing temperature and field. 

Observation of a supercurrent and multiple Andreev reflections in an SNS device requires 
a normal channel length L shorter than two principal length scales: the thermal length 

which characterizes the spatial decay of the pair amplitude in the N channel and the 
electron phase-breaking length L^. At low temperatures, the cut-off is typically determined 
by L^. For diffusive channels {L > /gi), ^n= ^/hD/2TTkBT, while for ballistic channels (L < 
^ei), = hvp/2'7ikBT. Here, D is the diffusion constant, vp is the Fermi velocity, ks is the 
Boltzmann constant and T is temperature. Since we do not directly know the carrier density 
n in our samples, we estimate lei and based upon some reasonable assumptions. In the 
normal state, the magnitudes of the resistivity and magnetoresistance of this nanoribbon 
(~ 0.51mf2-cm) are similar to measurements in Bi2Se3 thin filmd^ with a carrier density 
n ~ 10^^ cm^. We thus use this carrier density as an upper limit for n and use the Drude 
model (Zei = vprnes / pne"^) to yield le\ ~ 50 nm i.e. transport in the channel must be 
diffusive. Using the experimentally reported effective mass for Bi2Se3, rries ~ O.lSmg, at 
roughly this carrier density and for transport normal to the c-axi^, we estimate a Fermi 
velocity f p ~ 4 x 10^ ms~^ from vp = ^^{Sn'^nY^^ and a bulk diffusion constant D ~ 
7 X 10"^ m^s-i from D = vphi/^- This yields ~ 70 nm at T = 1.8 K. We note that 
even if the assumed nanoribbon carrier concentration was an order of magnitude lower (i.e. 
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10^^ cm~^), the corresponding mean free path lei would be ~ 250 nm and thermal length 
~ 100 nm at 1.8 K, both of which are still smaller than the channel length. Next, 
we estimate a lower bound for by extracting the electron inelastic length from weak 
antilocalization fits to the magnetoconductance of Bi2Se3 nanoribbons with normal metal 
(Pt) contacts. (The large He for the W contacts precludes a direct measurement of the low 
temperature magnetoconductance in samples with superconducting contacts.) In the two 
devices measured (Fig. 3c and 3d), we find 140 nm at T = 2 K. 

Thus, our analysis shows that we are observing a supercurrent and multiple Andreev 
reflections in SNS junctions with N channel lengths that are significantly longer than the 
values for and deduced from the diffusive normal state transport. The simplest expla- 
nation for this apparent inconsistency is the coexistence of ballistic (surface) and diffusive 
(bulk) transport channels. We speculate that the latter may dominate the transport in the 
normal state, but that the proximity-induced superconductivity preferentially occurs in the 
ballistic surface channel. The validity of this conjecture can be checked by estimating in 
the ballistic limit, using the value of the Fermi velocity f p ~ 5 x lO^ms"^ for the TI surface 
state in Bi2SeJ^ this yields a value ^N{ciean) = 2nkBT ~ at T = 1.8 K, which is more 

consistent with the observation of a supercurrent at comparable channel lengths. 

Our hypothesis of a ballistic surface channel that preferentially supports proximity- 
induced superconductivity is bolstered by transport measurements in three other 
Bi2Se3 nanoribbon devices. Figure 5 presents data taken on device B, a Bi2Se3 nanorib- 
bon of thickness t = 60 nm and width w = 430 nm. The normal state resistivity of this 
nanoribbon (1.37 mr2-cm) is almost three times greater than device A, yielding shorter 
estimated elastic and thermal lengths {lei ~ 20 nm, ~ 40 nm). This device has two mea- 
surable channels (Bl and B2 in the inset to Fig. 5a) with conservatively estimated channel 
lengths Lbi ~ 0.44 yum and Lb2 ~ 1.05 /im, respectively. Figures 5a-b show the onset of 
the proximity effect in these channels, while Figs. 5c-d show the bias-dependent differential 
conductance. The presence of multiple Andreev reflections in channel Bl and their absence 
in channel B2 are both consistent with our earlier estimate of ^N(ciean) ~ 340 nm for a clean 
surface channel since L^i ~ 'CN(cican) ^ Lb2- We have also observed subharmonic gap struc- 
ture in two additional nanoribbons (not shown). The overall behavior is similar to that 
discussed above. 
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V. TRANSPORT MEASUREMENTS WITH SUPERCONDUCTING (TUNG- 
STEN) CONTACTS: ANOMALOUS MAGNETORESISTANCE OSCILLATIONS 

We now discuss the occurrence of unusual MR oscillations in Bi2Se3 nanoribbons under 
conditions where the contacts become superconducting but where a supercurrent has not 
yet been established. In order to make consistent comparisons, we carried out both sets of 
measurements in the two probe geometry on device B, recalling that the electrode configu- 
ration shown in Fig. 5a precludes a true four probe measurement for the shorter channel. 
As a check, we also measured channel B2 in the four probe geometry and found identical 
behavior as the two probe measurements. Figure 6a shows the MR for channel Bl at T = 
2.2 K with a magnetic field applied perpendicular to the plane of the nanoribbon. At low 
field {\H\ < 0.6 kOe), we observe a sharp cusp characteristic of weak anti-localization, sim- 
ilar to earlier results from Bi2Se3nanoribbons contacted with normal electrodes. However, 
at higher fields, we find well defined and reproducible MR oscillations superimposed on a 
smooth background. The inset to Fig. 6a plots the magnetic field position of both peaks 
and valleys. If we exclude the low field region, we find that the MR oscillations are periodic 
in H with a period AH ~ 3.22 kOe. Figure 6b shows that the oscillation amplitude varies 
non-monotonically with temperature, reaching a maximum in the vicinity of T = 2.2 K and 
vanishing at both lower and higher temperatures. Similar behavior is observed in channel 
B2 except with a smaller period AH ~ 1.44 kOe (Figs. 6c and 6d). Similar oscillatory MR 
has also been observed in 2 other devices. For example. Fig. 7 shows the MR of device 
D in three primary field orientations: perpendicular to the plane of the nanoribbon (Fig. 
7a and 7b), parallel to the nanoribbon axis (Fig. 7c), and parallel to the plane but per- 
pendicular to the axis of the nanoribbon (Fig. 7d). While similar MR oscillations that are 
periodic {AH ~ 1.60 kOe) at higher field were observed when the field is perpendicular to 
the nanoribbon plane, they were not observed in the other two primary field directions. 

The observed MR oscillations cannot be explained using standard scenarios. The 
Shubnikov-de Haas effect is ruled out since the oscillations are not periodic in 1/H. The 
Aharonov-Bohm effect of surface currents-- is also ruled out for several reasons: the os- 
cillations are only observed when the magnetic field is perpendicular to the nanoribbon 
plane and not when it is parallel; the oscillation amplitude has a non-monotonic tempera- 
ture dependence; finally, MR oscillations are absent when using normal metal contacts (e.g. 
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Pt down to T = 2 K or W contacts above the superconducting transition temperature), 
clearly indicating that important role of superconductivity. We note that the Aharonov- 
Bohm oscillations were not observed in device D in parallel field due to the dominance of 
the normal transport in a long channel {L ^ ^N{ciean)) even when the electrodes are in the 
superconducting state. We also discount another conventional explanation for the oscilla- 
tory MR, namely fiuxoid quantization which requires a multiply connected geometry where 
fiux can penetrate normal regions encompassed by a superconductoi'^. Assuming fiux- 
oid quantization, the observed oscillation period implies a characteristic geometrical area 
A = {h/2e){AH)-^ = 0.643 x 10^ nm^ and 1.44 x 10^ nm^ for channels Bl and B2, respec- 
tively. This is far smaller than any relevant geometrical regions in our devices. The scaling 
of the oscillation period with the length of the Bi2Se3 channel (Bl and B2), roughly doubling 
when the channel length is halved, further rules out any insidious fiuxoid quantization that 
might be associated with FIB damage in the contact region. 

Having dismissed standard explanations for our observations, we resort to a new model 
that predicts MR oscillations when a thin, narrow superconducting strip is placed in a per- 
pendicular magnetic fielcP^'^. We are cautiously aware that this 'Weber blockade" model 
was developed for a genuine superconductor, while we are assuming that it can also be 
applied to a normal metal supporting proximity-induced superconductivity. There is the- 
oretical basis for the formation of vortex-like pair correlations in the N region of diffusive 
SNS junctional, but we are unaware of a similar prediction for a ballistic case such as the 
proximity effect in a high mobility TI surface state. In the Weber blockade model, the per- 
pendicular magnetic field controls the number of Pearl vortices in the superconducting strip. 
Maximum dissipation (i.e. maximum resistance) occurs when the energies of configurations 
with and A^ + 1 vortices are degenerate, and the MR then has a period related to the 
sequential addition of vortices with increasing magnetic field. Certain observations are im- 
mediately consistent with this model: for instance, the onset of periodic oscillations occurs 
only after a finite field (see insets to Figs. 6a and 6c) and the oscillation period is larger than 
that expected using fiuxoid quantization in the geometric area of the channel. Assuming 
that a proximity induced superconductor can be described by a uniform superconducting 
coherence length ^ (distinct from the length scale of the diffusion of Cooper pairs into the 
superconductor), this model can be extended to predict ^ in the region of proximity induced 
superconductivity for our geometry. 
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Using the geometrical parameters for our nanoribbons, we calculated the period of the 
MR oscillation as a function of the coherence length ^ (i.e. the size of a vortex). Figure 
6e shows the result of this calculation for channels Bl and B2, with channel width of 430 
nm and nominal channel lengths of 0.94 /im(squares) and 1.55 /im(circles), respectively. A 
comparison with the experimentally observed period yields ^ = 191(186) nm for channel 
Bl (B2). This length scale corresponds to the size of a vortex in the Bi2Se3 channel and, 
as expected, is much larger than the coherence length of a typical superconductor (such as 
the tungsten contacts). The inset to Fig. 6e shows, in a model calculation, how vortices 
distribute in a nanoribbon as the energy (magnetic field) increases. In Fig. 6f, we plot the 
number of vortices in the channel as a function of the applied magnetic field, assuming the 
fitted values of ^ in Fig. 6e. The model assumes a uniform superconducting gap; this may 
account for the disagreement between the experiment and the calculation using the fitted 
value of ^ for the first few vortices (Fig. 6f). 

VI. CONCLUSIONS 

In summary, we have demonstrated proximity-induced superconductivity in 
Bi2Se3 nanoribbons contacted with superconducting electrodes. Our experiments 
suggest that even when diffusive bulk conduction coexists with ballistic surface conduction 
in a 3D TI, the proximity effect may preferentially couple to the latter. Furthermore, 
the superconductor-TI device configuration demonstrated here provides a viable route for 
creating vortices near the interface between these two classes of materials. The formation 
of these vortices is manifest in MR oscillations whose period can be well explained using a 
simple model. Thus, our experiments show that superconductivity - and possibly vortices 
- can be realized in experimental geometries that are directly relevant to the search for 
Majorana fermions in condensed matter systemd^. 
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TABLE I: Summary of Bi2Se3 nanoribbon devices 



Device Length Width Thickness Contact width Resistivity at 6 K Contact material 
()Um) (nm) (nm) (nm) (mO-cm) 



A 


1.08 


600 


60 


200 


0.51 


W 


Bl 


0.94 


430 


60 


430 


1.37 


w 


B2 


1.55 


430 


60 


430 


1.37 


w 


C 


2.47 


284 


276 


316 


0.79 


w 


D 


2.30 


240 


77 


240 


0.70 


w 


E 


5.12 


300 


98 


122 


0.69 


w 


F 


1.00 


700 


60 


300 


3.53 


w 


G 


5.73 


270 


60 


234 


0.39 


Pt 


H 


0.94 


370 


50 


280 


0.48 


Pt 
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Raman shift (cm ) 

FIG. 1: (Color online) a. SEM image of as-grown Bi2Se3 nanoribbons on a Si substrate, b. 
HRTEM image from the edge of a Bi2Se3 nanoribbon. The interplanar distance along the growth 
direction [1120] is 0.21 nm. c. Selected area electron diffraction pattern from the nanoribbon in 
Fig. lb shows hexagonal symmetry. The growth direction is along [1120]. d. A typical Raman 
spectrum from a single Bi2Se3 nanoribbon with 514.5 nm excitation at room temperature. All 
three vibration modes are consistent with those from bulk Bi2Se3. The inset shows an optical image 
of a nanoribbon supported over a hole of a TEM grid. The green spot is the laser illumination. 



15 



500 1000 

X (nm) 



FIG. 2: (Color online) a. High-angle annular dark field image of a W strip deposited on a 
Bi2Se3 nanoribbon supported on a holey carbon TEM grid using the same deposition conditions 
for the actual devices. In the region next to the W contact (indicated by the circle), selected area 
diffraction pattern (inset) shows clear hexagonal Bi2Se3 single crystal structure. The line on the 
nanoribbon indicates the position of the X-ray energy dispersive line profile and the arrow shows 
the line scan direction. EDS spectra separated by 50 nm were collected in STEM mode using drift 
correction, b. EDS counts (arb. units) vs. position plot of elemental W, Bi, and Se. The W signal 
is within the background at ~ 450 nm, which corresponds to a spread of ~ 250 nm from the left 
visible edge of the tungsten contact. 
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FIG. 3: (Color online) a. Temperature dependent conductivity for device G. The solid line is 
a ln(T) fit. b. Temperature dependent conductivity for device H. The solid line is a ln(T) fit. 

c. Magnetoconductivity for device G at T = 3.0 K (squares) and T = 2.0 K (circles). The solid 
lines are fits to the Hikami-Larkin-Nagaoka theory with electron-electron interaction included. 

d. Magnetoconductivity for device H at T = 3.0 K (squares) and T = 2.0 K (circles). The solid 
lines are fits to the Hikami-Larkin-Nagaoka theory with electron-electron interaction. 
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FIG. 4: (Color online) a. Two probe zero bias dV/dl vs. temperature for device A at = 0. 
Inset shows an SEM image of the device, with the arrow indicating the measured channel with 
edge-to-edge length of 1.08 ^m. between two W electrodes, b. I — V characteristics of device A 
at various temperatures, measured using the same contacts as in (a), c. dl/dV vs. V in device 
A at T = 500 mK and in zero magnetic field. The arrows identify a consistent subharmonic series 
of conductance anomalies corresponding to subharmonic gap structure with n = 2,4, 8). d 
Position of differential conductance anomalies as a function of the index 1/n. e Surface plot of 
dl/dV as a function of the current / and temperature at H = 0. f Surface plot of dl/dV as a 
function of the current / and perpendicular magnetic field H at T = 500 mK. 
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FIG. 5: (Color online) a. Two probe resistance vs. temperature for channel Bl in different 
perpendicular fields. Inset shows an SEM image of the device B, with the two arrows indicating 
channels Bl and B2 with edge-to-edge length of 0.94 ^um and 1.55 fim, respectively, between two 
W electrodes, b. Four probe resistance vs. temperature for channel B2 in different perpendicular 
fields, c. dl/dV vs. V for channel Bl at different temperatures and in zero magnetic field using 
a two probe configuration. Data for different temperatures are shifted or scaled with respected 
to the curve at T = 4.0 K for clarity. The arrows indicate conductance peaks corresponding to 
subharmonic gap structure (2A/ne with n = 1,2,6). d. dl/dV vs. V for channel B2 at different 
temperatures and in zero magnetic field using a two probe configuration. The data for different 
temperatures are shifted or scaled with respected to the curve at T = 4.0 K for clarity. 
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FIG. 6: (Color online) a. MR in device B (channel Bl, see lower inset) at T = 2.2 K. The upper 
inset plots the position of both peaks and valleys, b. MR in channel Bl at different temperatures. 
The data are shifted with respect to the data at T = 2.2 K for clarity, c. MR in channel B2 
at r = 2.6 K. The inset plots the position of both peaks and valleys, d. MR in channel B2 at 
different temperatures. The magnetic field in panels a-d is perpendicular to the nanoribbon plane, 
e. Calculated MR oscillation period vs. coherence length (size of vortex) for channels Bl (squares) 
and B2 (circles). Comparison with the data in Figs. 2a and 2c yields ^ = 191 nm and ^ = 186 
nm for channel Bl and B2, respectively. The inset is a model calculation showing how vortices 
distribute in a nanoribbon as the energy (magnetic field) increases, f. Magnetic field vs. number 
of vortices. Open symbols are experimental data and filled symbols are fits using the values of ^ 
obtained in panel e. 
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FIG. 7: (Color online) a. MR for device D at T = 2.2 K in field perpendicular to nanoribbon 
plane using four probe geometry. Pronounced and reproducible MR oscillations with a period of 
AH ~ 1.6 kOe for \H\ > 5.5 kOe were observed. The inset is the FFT of the MR after a smooth 
background is subtracted for H > 5.5 kOe. b. MR at different temperatures for device D in field 
perpendicular to the plane of the nanoribbon. c. MR at different temperatures for device D in 
field parallel to the axis of the nanoribbon. d. MR at different temperatures for device D in field 
parallel to the plane but perpendicular to the axis of the nanoribbon. 



21 



